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ABSTRACT 

A  coupled  air-ice-ocean  theoretical  model  for  the  marginal  ice  zone  (MIZ)  is  developed  and  used  to  show 
that  an  off-ice  and  divergent  wind  field  not  only  produces  a  dilation  of  the  MIZ  (as  is  generally  thought),  but 
also  generates  a  compaction  of  MIZ  for  some  circumstances.  An  ice  divergence-convergence  criterion  in  the 
MIZ  is  found  and  used  to  help  explain  some  features,  such  as  ice  edge  bands,  and  the  formation  and  maintenance 
of  polynyas  and  leads. 

The  model  contains  three  parts:  a  thermally  forced  boundary  layer  air  flow,  a  mechanically  driven  ice  drift, 
and  a  reduced  gravity  ocean.  The  three  components  are  linked  through  the  surface  temperature  gradient  and 
various  interfacial  stresses.  Model  results  show  that  in  the  MIZ,  a  thermally  generated  surface  wind  blowing 
from  ice  to  water  (ice  breeze)  is  generally  divergent  over  ice  owing  to  the  maximum  surface  temperature  gradient 
located  at  the  ice  edge.  Such  divergent  local  winds  can  make  the  ice  flow  either  divergent  or  convergent,  depending 
on  the  properties  of  ice  and  water  underneath.  For  thin  ice  and  a  relatively  deep  surface  water  layer,  the  ice 
turning  angle  (angle  between  the  direction  of  surface  wind  and  ice  flow)  is  small  and  the  ice  flow  is  divergent. 
In  contrast,  for  thick  ice  and  a  relatively  shallow  surface  water  layer,  the  turning  angle  is  large  and  the  ice  flow 
is  convergent. 


l.  Introduction 

In  this  study,  a  steady  state  coupled  air-ice-ocean 
model  is  employed  to  investigate  an  ice  breeze  mech¬ 
anism  for  ice  divergence/convergence  in  the  MIZ.  The 
importance  of  ice  breeze  in  the  MIZ  is  shown  as  follows. 
Schmidt  (1947)  used  a  simple  linear  model  to  discuss 
seabreeze/landbreeze  phenomena  and  concluded  that 
for  a  maximum  land-water  temperature  gradient  of 
4.3°C/100  km,  a  maximum  landbreeze  intensity  of 
about  2ms-1  would  be  reached.  From  observations 
in  the  Southern  Bering  Sea,  Reynolds  et  al.  (1985)  es¬ 
timated  that  ice  floes  drift  to  the  right  of  the  wind  by 
approximately  30° 'at  about  4%  of  the  wind  speed  at  3 

m.  The  ice  breeze  is  analogous  to  the  landbreeze. 
Therefore,  1°C/100  km  surface  temperature  gradient 
roughly  produces  0.5  m  s-!  ice  breeze,  which  in  turn 
generates  about  2  cm  s'1  ice  flow.  The  CTD  data  in 
the  Central  Bering  Sea  ice  edge  (Muench,  1983)  show 
that  the  cross  ice  edge  sea  surface  temperature  (SST) 
gradient  is  around  7°C/100  km.  However,  the  mete¬ 
orological  observations  in  the  same  area  (Lindsay  and 
Comiskey,  1981)  indicate  that  the  surface  air  temper¬ 
ature  gradient  across  the  ice  edge  is  nearly  5°-10°C/ 
100  km.  Therefore,  near  the  Bering  Sea  ice  edge,  the 
ice  breeze  is  around  2.5-5  m  s-1  (Schmidt’s  estimation), 
and  the  associated  ice  floe  drift  is  about  10-20  cm  s-1 
(Reynold’s  estimation). 

During  September-October  1979  the  Norwegian 
Remote  Sensing  Experiment  was  carried  out  in  the 


MIZ  north  of  Svalbard.  It  was  found  that  near  ice  edge 
the  wind  generally  blew  from  ice  to  ocean  (except  dur¬ 
ing  26-27  September)  and  the  ice  drift  vector  was  de¬ 
flected  to  the  right  of  the  wind,  with  a  further  deflection 
of  the  current  in  the  mixed  layer  (Johannessen  et  al., 
1983). 

A  possible  mechanism  for  the  mesoscale  air-ice- 
ocean  interaction  in  the  MIZ  is  depicted  in  Fig.  1 .  The 
low-level  air  flow  generated  by  a  surface  temperature 
gradient  moves  from  ice  to  water  with  some  deflection 
to  the  right  (or  left)  in  the  north  (or  south)  polar  region. 
This  ice  breeze  should  have  a  maximum  speed  near 
the  ice  edge  and  decrease  both  iceward  and  waterward, 

1. e.,  it  is  divergent  over  the  ice  and  convergent  over  the 
water.  The  ice  flow  driven  by  the  surface  wind  field  is 
to  the  right  (or  left)  of  the  wind  in  the  north  (or  south) 
polar  region.  The  turning  angle  between  ice  flow  and 
surface  air  flow  varies  from  0°  to  90°,  depending  on 
the  properties  of  the  ice  and  water  underneath. 

The  air-ice-ocean  interaction  model  depicted  in  the 
subsequent  sections  is  used  to  determine  the  main  cri¬ 
terion  of  ice  divergence/convergence. 

2.  Thermally  forced  boundary  layer  air  flow 

In  this  section  we  utilize  the  same  planetary  bound¬ 
ary  layer  model  treated  by  Kuo  (1973)  and  Chu  (1986, 
1987a,b,c,d)  to  simulate  a  thermally  forced  boundary 
layer  air  flow.  The  coordinate  system  is  chosen  such 
that  the  x-axis  is  in  the  cross  ice  edge  direction  and  the 
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where 

£/=g8D7y(2LO0*)  (4) 

is  the  scale  of  icebreeze  wind.  If  we  assume  that  the 
local  air  flow  satisfies  the  modified  Boussinesq  approx¬ 
imation  (Kuo,  1973),  the  vorticity  equation,  the  mo¬ 
mentum  equation  (both  in  y  direction),  and  the  heat 
equation  for  the  air  disturbances  generated  by  differ¬ 
ential  surface  temperature  gradient  near  the  MIZ  are 
(Chu,  1986,  1987a,b,c,d) 


Water  Current 


Surface  Temperature  Gradient 
Cold  Ice  and  Warm  Water 


I 


Icebreeze  Circulation 


' — [ice  Flow 


Water  Current 


Fig.  1 .  Ice-air-ocean  coupled  system. 


y-axis  parallels  the  ice  edge,  as  shown  in  Fig.  2.  The  x- 
coordinate  unit  length  is  twice  the  MIZ  width,  L  =  200 
km,  and  that  of  the  vertical  coordinate  is  5  =  (e/0)i/2, 
where  v  is  the  vertical  eddy  viscosity,  and  ft  the  angular 
velocity  of  the  earth’s  rotation.  The  line  x  =  'h  is  located 
at  the  ice  edge.  It  is  considered  that  spatial  variations 
in  the  MIZ  are  much  larger  perpendicular  to  the  ice 
edge  than  parallel  to  it,  and  hence  derivatives  with  re¬ 
spect  to  y  are  assumed  to  be  zero  (Lepparanta  and  Hib- 
ler,  1985). 

The  potential  temperature  of  air  is  divided  into  two 
parts:  a  basic  state  0b*(z)  and  perturbation  z).  The 
basic  state  is  given  by 

0e*(z) =  0bo»  +  (N26m  m&/g)z,  ( 1 ) 

where  the  subscript  means  dimensional  quantities, 
0B o*  is  the  basic  air  potential  temperature  at  the  surface, 
and  AThe  Brunt- Vaisala  frequency.  Observations  show 
that  the  surface  air  potential  temperature  perturbation 
0*(x,  0)  monotonically  decreases  from  water  to  ice.  We 
use  the  following  sinusoidal  function  to  parameterize 
the  surface  thermal  condition: 

ff*(x,  0)  =  —  DT0  cosirjt  JL,  (2) 

where  DTo  is  the  characteristic  temperature  difference 
across  the  MIZ.  The  coordinates  and  atmospheric 
variables  are  nondimensionalized  by  setting 

(a*  ,  z*)  =  (xL,  zS)  ) 


s'  =  6'*/0m  =  (DTo/6m)S^ 

(u(? ,  »<“> ,  w^)  =  [/(M(a),  via\  w(a)o/L)  > 
P(?  =  (g8DT0/8m)p(a) 

s  =  (p/ny/2, 


(3) 


Ed*&a)/dz4  =  fodv(a)/dz  - 

ds<a)/dx. 

(5) 

Ed2v{a)/dz2  =  -fod^/dz. 

(6) 

Ed2s(a)/dz2  =  Rid&a)/dx. 

(7) 

In  the  foregoing, 

u(a)  = -dV(a)/dz,  w(a)  =  d^/dx. 

(8) 

where  f0  =  simp,  <p  is  the  latitude,  and 

£=  v/(2Q52)  =  ‘/zj 
Ri  =  (5N/2LQ)2  J 

(9) 

are  the  Ekman  and  Richardson  numbers,  respectively. 
Eliminating  v(a)  and  sXa)  from  (5)-(7)  we  find  that  the 
streamfunction  satisfies  the  following  partial  differential 
equation: 

d  4/dz4  +/02Jd2*(a)/dz2  +  R\d2¥a)/dx2  =  0.  (10) 

We  solve  ( 1 0)  for  the  streamfunction  'R<"),  and  obtain 
the  solutions  of  v(a)  and  s(a)  from  (6)  and  (7)  after  sub¬ 
stituting  \k(fl).  The  local  air  flow  is  thermally  forced  by 
the  surface  temperature  gradient  as  indicated  in  (2), 
therefore  the  streamfunction  should  be  written  as 

^(a)(A,  z)  =  *la)(z)  sinTTX.  (11) 

The  boundary  conditions  in  the  vertical  direction  are 
derived  as  follows.  The  dependent  variables  should  re¬ 
main  finite  as  z  -*■  oo,  i.e., 
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lim(|tf(fl)|,  \d&a)ldz\,  |u(fl)|,  |i(a)|)  <oo.  (12) 

z-»oo 

The  boundary  conditions  at  z  =  0  are 

^(a)  =  0,  d^aVaz  =  Md2-¥a)ldz2, 

v(a)  =  Mdv(a)/dz,  sia)  =  — coswx,  ( 1 3) 

where  M  is  a  measure  of  the  effective  depth  of  the  con¬ 
stant  stress-sublayer  (Kuo,  1973).  Substituting  (1 1)  into 
(10)  we  obtain  the  following  sixth  order  ordinary  dif¬ 
ferential  equation  for  the  Fourier  coefficient  Sk<a): 

d6*{a)/dz6  +  4tfd2va)ldz2  -  4ir2  Ri^(a)  =  0.  (14) 

The  general  solution  of  (14)  has  the  following  form: 

6 

*(a)(z)  =  2  aj  exp(Xyz),  (15) 

7=1 

where  the  eigenvalues  Xj(j  =  1,2,  •  •  •,  6)  are  the  roots 
of  the  sixth-order  algebraic  equation: 

X6  +  4/o2X2-47T2Ri  =  0.  (16) 

According  to  the  upper  boundary  conditions  listed  in 
(12)  we  must  set  up  coefficients  which  correspond  to 
those  eigenvalues  with  positive  real  parts  equal  to  zero. 
Consequently  the  general  solution  (16)  satisfying  the 
top  boundary  conditions  may  be  written: 

*ia\z)  =  2  tfy  exp(XjZ),  (17) 

7=1 

where  the  eigenvalues  X,  all  have  negative  real  parts. 
Substituting  (17)  into  ( 1 1 )  we  obtain  the  streamfunction 

3 

'k(a)(x,  z,  t)  =  2  aj  exp(Xyz)  simrx.  (18) 

7=  i 

Integrating  the  momentum  equation  (6)  and  the  heat 
equation  (7)  with  respect  to  z  after  substituting  (18), 
we  find  that  v(a)  and  s(a)  are  given  by 

3 

v(a)(x,  z)  =  b-2/0Z ( aj/Xj )  exp (X,z)  simrx,  ( 1 9) 

7=1 

3 

s(a)(x,  z)  =  2x  Ri  2  (£>/X/)  exp(Xyz)  cos7rx.  (20) 

7=1 

Substituting  the  solutions  (18)-(20)  into  the  surface 
boundary  conditions  listed  in  (13)  we  obtain  the  fol¬ 
lowing  four  algebraic  equations  for  a,  and  b: 

3 

2«,  =  0,  (21) 


'ZXJ(l-M\j)aj  =  0,  (22) 

7=1 

3 

-2/o  2  (fly/XyX  1  -  M\j)  +  6  =  0,  (23) 

7=1 


3 

2t  Ri  2  <7/ A/  =  —  1 .  (24) 

7=1 

The  four  constants  at ,  a2,  a3  and  b  can  be  readily  ob¬ 
tained  by  solving  the  four  linear  nonhomogeneous  al¬ 
gebraic  equations. 

The  surface  wind  (dimensional  form)  driven  by  the 
horizontal  temperature  gradient  is  computed  by 

u(£Xx,  0)  =  M0* (a)  simrx,  (25) 

v^(x,  0)  =  r  o* <a)  sinxx,  (26) 

where 

3 

m0.  (a)  =  u  Re(  2  a,  Xj)  (27) 

7=i 

3 

v0.<a> -  U Re[6  -  2/0  2  (aj/Xj)].  (28) 

7=1 

The  purpose  of  using  the  nondimensional  form  for 
the  atmosphere  is  to  get  solutions  that  simply  show  the 
magnitude  of  the  boundary  layer  air  flow.  After  ob¬ 
taining  solutions  (27)  and  (28),  we  will  use  dimensional 
variables.  The  subscript  will  be  dropped.  Figure  3 
shows  the  cross-ice-edge  atmospheric  circulation 
(u(a\  \da))  at  <p  =  75°N,  DT0  =  6  K  (mean  temperature 
gradient  ttDTq/L  =  9.4°C/100  km)  and  N  =  1(T2  s“‘. 
We  see  the  obvious  ice  breeze  at  low  levels.  The  speed 
of  the  ice  breeze  reaches  a  maximum  (which  is  around 
4  m  s-1)  at  the  ice  edge,  which  is  consistent  with 
Schmidt’s  (1947)  results. 

3.  Free  ice  drift  model 

Roed  and  O’Brien  (1983)  show  that  the  internal  ice 
stress  and  the  nonlinear  terms  are  not  crucial  for  the 
MIZ  dynamics,  and  that  the  ice  velocity  is  an  order  of 
magnitude  larger  than  water  current  in  the  MIZ.  Fol¬ 
lowing  their  theoretical  model,  the  momentum  equa¬ 
tions  for  a  free  drift  ice  model  are 

fv{,)  +  [txm  +  Tx(^]/(p,//,)  =  0,  (29) 

-fu(i)  +  [ryM  +  =  0,  (30) 

where  /  (=  2fi  sirup)  is  the  Coriolis  parameter,  V(,) 
=  (m(,),  t>(,))  the  ice  velocity,  H,  the  mean  ice  thickness, 
Pi  the  ice  density,  T(a,,)  =  the  wind  stress 

on  ice,  r(K’,')  =  (tx(w’’\  Ty(Wf,))  the  water  stress  on  ice. 
The  air  and  water  stresses  on  ice  are  computed  by 

riaJ)  =  PaCat V(fl)  (31) 

=  pwCwi[V(M’>-V(')],  (32) 

where  Cai  and  Cwi  are  dimensional  (ms-1)  air  and  water 
drag  coefficients  (on  ice);  pa  and  pw  are  the 
densities  of  air  and  surface  water,  respectively;  V(a) 
=  (u<a\  v(a))  is  the  surface  wind,  and  V(M,)  =  (u(wi, 
v(W))  the  water  velocity. 
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Fig.  3.  Thermally  forced  ice  breeze  circulation  in  the  cross-ice 
section  with  u  in  m  s-1  and  w  in  cm  s-1. 

4.  Reduced  gravity  ocean  model 

Suppose  that  the  ocean  is  composed  of  two  layers, 
in  which  the  horizontal  inhomogeneity  of  density  inside 
the  same  layer  is  much  smaller  than  the  density  dif¬ 
ference  between  The  two  layers,  and  the  lower  layer  is 
deep  enough  for  motion  in  the  lower  layer  to  be  van¬ 
ishingly  small.  Such  a  model  is  referred  to  as  a  reduced- 
gravity  model.  Furthermore,  the  model  has  the  simplest 
form  of  dissipation,  namely,  Rayleigh  friction  with  de¬ 
cay  rate  R  and  Newtonian  cooling,  also  with  decay  rate 
R.  The  momentum  and  continuity  equations  for  the 
ocean  are  written  by  (Gill,  1982;  Roed  and  O’Brien, 
1983) 

Ru(w)—fv(w)  =  -g*dh{w)/dx 

+  [( 1  -  A)tx + Arx^]/(pwHw)  (33a) 

Rv(w)  +fu(w)  =  [(1  —  A)ry(a-W)  +  AT«’w)y(PwHw)  (33b) 

Rhiw)  +  Hwdu{w)/dx  =  0,  (33c) 

where  g*  is  a  reduced  gravitational  acceleration,  Hw 
the  equilibrium  depth  of  the  pycnocline,  h{w)  the  thick¬ 
ness  of  the  upper  layer,  and  =  (rx(,,K,),  Tyl,’w))  the 
ice  stress  on  the  water.  has  the  same  magnitude 
but  the  opposite  direction  as  t(wj\  i.e., 

?('»  =  -YwJ)  =  -pwCwl(Yw)  -  V(/>).  (33d) 

Here  A  is  the  ice  compactness  or  ice  concentration 
(i.e.,  the  fraction  of  area  covered  by  ice).  Hibler  and 
Ackley  (1983)  took  the  50%  limit  of  A  as  the  ice  edge. 
Inside  the  MIZ,  the  ice  concentration,  A,  decreases  wa- 
terward  very  slowly  from  the  boundary  between  pack 
ice  and  the  MIZ  (where  A  ~  1)  to  some  place  near  the 
ice  edge,  and  then  diminishes  very  rapidly  to  the  ice 
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Table  1 .  The  standard  model  parameters. 


L  =  200  km 

»  =  Sm!s"' 

»*=  10s  m2s-' 

SI  =  0.7292  X10-4s-‘ 

g  =  9.81  ms-2 

R  =  2.08X  10— 5  s~ 1 

e«,>  =  270K 

CWI  =  0.86X 

Cat  =  0.03  m  s‘l 

pa  =  1 .29  kg  m"3 

pw=  1000  kg  m-3 

Pi  =  9 1 0  kg  rrT3 

g*  =  0.005  m  s-2 

D7o  =  6K 

edge.  If  the  area  we  focus  on  is  not  very  close  to  the 
ice  edge,  we  may  set  A  ~  1 .  The  portion  of  air  stress 
on  the  water,  ( 1  -  A)T(a’w),  is  then  neglected  compared 
to  the  ice  stress  on  the  water,  At(,'w).  Setting  A  =  1  and 
eliminating  hM  from  (33a)  and  (33c),  the  momentum 
equations  (33a)  and  (33b)  become 

Ruiw)  -  fv{w)  -  (g*H  JRL2)d2u(w)/dx2  =  tx(,'w)/(pwHw) 

(34) 

Rv(K)  +MW)  =  ry^/(pwHw).  (35) 

5.  Ice  divergence-convergence  criterion 

The  surface  wind  is  the  forcing  term  of  the  ice-ocean 
system.  Since  our  model  is  linear,  the  mechanically 
forced  ice  and  water  motion  V(<>  and  V(H,)  should  have 
the  same  Fourier  component  as  the  forcing  term  V(a), 
ie., 

u(1)  =  UoU)  sin7rx,  u(,)  =  i>0(,)  sinirjc,  (36) 

u(w)  =  u0iw)  sin-irx,  v(w)  =  v0(w)  sin7rY.  (37) 

Substituting  (31),  (32),  (33d),  (36)  and  (37)  into  (29), 
(30),  (34)  and  (35)  we  have  the  following  four  linear 
algebraic  equations  for  the  Fourier  coefficients  uQ(w\ 
v0{w\  UoU)  and  v0U): 

[R  +  Tr2g*/(RL2)  +  a>3Ww)  -fv0{w)  -  aW°  =  0,  (38) 
+  (R  +  -  oW°  =  0,  (39) 

-co2Mo(w)  +  co2Mo(i)  ~fv  0(i)  =  wlUoia)  (40) 

-oW1*0  +  co2Uo(i)  +fuo{0  =  o>iVoia\  (41) 


Rt 

Fig.  4.  Dependence  of  a  on  the  atmospheric  Richardson  number 
Ri  for  three  different  latitudes:  (a)  65°,  (b)  70°  and  (c)  75°. 
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Table  1  lists  the  parameter  values  for  the  present 
calculations.  Figure  4  describes  the  dependence  of 
the  ice  breeze  deflection  angle  (angle  between  sur¬ 
face  temperature  gradient  and  surface  wind)  a 
=  tan“1(Uo(a)/«o(‘J))  on  Ri  at  three  different  latitudes 
(65°,  70°  and  75°).  The  |a|  increases  with  increasing 
latitude  and  decreases  with  increasing  Ri.  At  65°  (N 
and  S)  latitude  it  varies  from  41.8°  when  Ri  =  10-3  to 
a  =  27.5°  when  Ri  =  4  X  1(T2. 

The  angle  between  surface  wind  and  ice  velocity,  (3 
(as  shown  in  Fig.  5),  is  defined  by 


where  parameters  a>(,  <o2  and  u3  are  defined  by 

Wi  =  PaCJiPjHj),  a>2  =  PwCJip.Hj),  co3  =  CJHW. 

(42) 

Solving  the  nonhomogeneous  linear  algebraic  equa¬ 
tions  (38)— (41)  we  obtain 

Mo(,)  =  «i[o>2(l  -  ot2)Uo(a)  +/(1  +  a3u2)v0{a)]/ 

[«22(1  ~  «i)(l  “  a2)+/(l  +  a3«2)2],  (43) 
U0(I)  =  aj,[<u2(l  -«iWa)-/(l  +a3 o>i)UoW]l 

[a)22(l  -  «,X  1  -  «2)  +/(1  +  «3^2)2],  (44) 
Mo(w)  =  +  w3)+/V°]/[«o(7?  +  «s)  +/2]  (45) 

Uo(w)  =  co3[a0Uo(,) -fu0(l)]/[a0(R  +  «3)  +/2]  (46) 

where 

a0^R  +  w3  +  Tc2g*Hw/(RL2) 

«i  =  co3(R  +  cu3)/[a0(R  +  w3)  +/2] 
a2  =  a)3a0/[a0(R  +  oo3)  +/2] 

«3  =  w3/[a0(R  +  w3)+/2].  (47) 


d  =  tan  l(v0M/uo(,))  —  a  (48) 

where  Wo(,)  and  v00>  are  computed  by  (43)  and  (44).  An 
ice  divergence-convergence  criterion  is  simply 

\a+J>ir12,  ice  convergence  (49) 

1<7t/2,  ice  divergence. 

We  compute  the  ice  turning  angle  /3  as  a  function 
of  Ht  and  Hw  for  <p  -  65°,  where  Ht  varies  from  0  to 
10  m  and  Hw  from  0  to  50  m.  The  distribution  of  0 
with  Hi  and  Hw  is  shown  in  Fig.  6.  Small  values  of  |/3| 
(small  turning  angles)  are  located  in  the  region  of  large 
Hw  and  small  //,  .  The  smaller  the  Ht  (or  the  larger  the 
Hw),  the  smaller  the  |/3|.  It  is  quite  understandable  that 
the  water  current  V<B,)  is  very  weak  as  Hw  becomes  large, 
due  to  small  mechanical  forcing  terms  &>3m(0  and 
co3u(')  in  (38)  and  (39).  When  the  water  current  is  weak 
(due  to  large  Hw)  the  turning  angle  largely  depends  on 
the  ice  thickness.  For  thin  ice,  Wi  P  f,  (40)  and  (41) 
show  that  V<0  deflects  by  small  angle  from  V(fl),  there¬ 
fore  the  ice  flow  is  divergent. 

Solving  the  algebraic  equation  (48)  for  Hi  after  setting 
|  a  +  /3|  =  7t/2,  we  get  the  dependence  of  Ht  on  Hw  for 


Hw  (M  ) 


FIG.  6.  Isolines  of  /3  (°)  in  the  (Ht,  Hw)  plane  for  <p  =  65°. 
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FIG.  7.  Ht-Hw  curves  for  zero  ice  divergence-convergence 
for  five  different  values  of  Ri. 


the  no  convergence  or  divergence  case,  as  shown  in 
Fig.  7.  The  curves  separate  the  parameter  plane  (//,, 
Hw)  into  ice  convergence  and  ice  divergence  parts.  We 
can  differentiate  between  ice  divergence  and  ice  con¬ 
vergence  due  to  the  ice  breeze  in  the  MIZ  by  the  pa¬ 
rameters  Ri,  Hi,  and  Hw.  If  the  drag  coefficient  Cwi  is 
doubled,  the  turning  angle  (3  only  has  minor  changes 
(less  than  10%). 

6.  Conclusion 

(i)  This  air-ice-ocean  model  is  intended  to  depict 
only  the  mesoscale  processes  of  air-ice-ocean  inter¬ 
action  in  the  MIZ.  The  synoptic-scale  pressure  gradient 
may  additionally  produce  surface  winds,  and  in  turn 
drive  ice  drift.  These  processes  are,  however,  beyond 
the  scope  of  this  paper.  Nevertheless,  when  the  ice-to- 
open-water  temperature  gradient  is  strong,  the  ice 
breeze  mechanism  discussed  here  may  become  as 
strong  as,  or  stronger  than,  the  synoptic  forcing. 

(ii)  The  model  shows  the  different  effects  of  ice 
breeze  on  ice  flow,  and  gives  an  ice  divergence-con¬ 
vergence  criterion  in  the  MIZ,  which  may  help  to  un¬ 
derstand  the  formation  of  ice  edge  bands,  and  the  for¬ 
mation  and  maintenance  of  polynyas  and  leads. 

(iii)  It  should  be  stressed,  however,  that  the  model 
cannot  be  used  to  describe  the  system  with  an  oceanic 
density  front  paralleling  the  ice  edge  as  observed  by 
Muench  and  Schumacher  (1985). 
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